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Rozwiązanie równania Fokkera-Plancka do rozszczepiajacych się gorących jąder 


1. INTRODUCTION 


The description of the fission process by a transport equation goes back 
to a fundamental paper by K. K ra mers [1] and has been taken up again by 
several authors (see e.g. [2-4]) in recent time, especially when one discovered 
that the transport theory was a useful tool to describe heavy-ions reactions. 

It was shown in [3] that within the linear response theory the fission 
dynamics is governed by a generalized Fokker-Planck equation (FPE) for 
the distribution function f(q,p,t): 


Of(a pt) _ josh. Of(q p,t) 8U(q) OF(a, pt) 
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Here m and y are the inertia and friction parameters respectively. U is the 
deformation potential which can be identified with the free energy in an 
adiabatic theory. 

The FPE can be solved analytically in the simplest cases only. It is 
the case for example when the transport coefficient is constant and the 
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potential U has a parabolic form. More complicated cases have to be solved 
numerically. We present in the next section some most frequently used 
numerical method of solving FPE. We shall also try to combine the Monte- 
-Carlo random numbers procedure with the propagator method. In the next 
section we present the numerical results for the probability that system does 
not fission. 


2. METHOD OF SOLVING OF FPE 


The approximate solution of the FPE with the variable transport 
coefficient for the fission dynamics was presented in [5]. It was assumed 
there that the form of the distribution function f(q,p,t) can be described 
by its first and second moments: 


1 
fla p, t) = Zadat $* 


x exp x OF -25p(- ON P) - Sy PY 


where Q and P are the average coordinate and momentum respectively: 


(2) 


P 


Q(t) = J qaf (4q, p, t)dqdp (3) 


P(t) = / Pf(q,p,t)dqdp (4) 


and $ is the matrix of the quadratic deviation from the average values Q 
and P: 


Sq = Jo - Q)” f (q, p,t)dgdp (5) 
Sæ = fla- Qp- P)4(a,,t)dadp (6) 
Spp = J (p - P)’ f (q, p,t)dqdp (7) 


The expression (2) is accurate for the gaussian form of the distribution 
function. Inserting eq. (2) into the FPE (1) one can get the coupled set of 
the differential equations [6] for the first: 


(8) 
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and the second moments: 
“aft = nów 00) 
os = Sy E Su - Hs, (11) 
z =2D- oo Se - LS, (12) 


This method is very convenient especially for the multidimensional case. 
The generalization of eqs (8-12) for that case is obvious [6]. But on the 
other side this gaussian approximation of the distribution function cannot be 
applied in the case of a large potential anharmonicity and variable transport 
coefficients. 

A better, more accurate method of solving of the FPE was proposed in 
[7] where the propagator method was applied. Namely it was assumed there 
that the density probability after a small time interval At is done by the 
folding integral from the density at the time ¢ and the propagator function 
K: 


f(a,p,t + At) = JEZ + At;q,p',t) f(q',p', t)dqdp (13) 


The propagator K was assumed to be of the gaussian form: 


K(q,p,t + At;q,p,t) = 


exp{— S229 1) = Spa = PP P')— SoolP— PY (44) 
2x det S : det $ 
where its widths S44, Sgp and Spp are equal to zero at the initial time t (it 
means that K is equal to the Dirac 6-function) and they grow with increasing 
time according to the same equations as those for the moments of the whole 
distribution, eqs (8-12). 

One has to repeat the process of evaluating of the folding integral (13) 
n times in order to get the final distribution at the time tn = to + nAt. 
The propagator method is accurate when the time step At is small enough. 
The numerical effort to proceed this recursive folding is very large and 
the calculation can be effectively done on large computers only. The 
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problem becomes even more complicated when we are dealing with the 
multidimensional collective space. This folding procedure becomes more 
difficult when one assumes the Gauss form of the distribution function in 
the moment direction. The validity of this approximation was discussed 
in [8]. 

In the present paper we propose a simplified version of the propagator 
method which nevertheless leads to the exact solution of FPE. Following the 
idea of [9] we adopt the Monte-Carlo method of generating random points 
on (q,p) plane to evaluate the folding integral (13). 


ty tyetorat tosty*Zat th? tornat 


initial dist. final dist. 


q 
Fig. 1. The schematic illustration of the Monte-Carlo method of solving of FPE 


The method of the calculation is illustrated in the graph in Fig. 1. and 
it is following: 

1. At the initial time tọ = 0 one has to generate N random points: 
(Po), where v = 1,2,...N. The distribution of these points is obtained 
from the initial shape of the function f(q,p, to) = 0. 

2. The position (qi;,py) of each random point and its smearing width 
S% (tı) at the time ty = to + At is given by the solution of the coupled 
differential equations (8-12). The initiał conditions for the widths are: 
Sqq(to) = Sgp(to) = Szp(to) =0 

3. In the next step, at the time t = t; only one random point (43, p3) 
for each trajectory v is generated. The random number generator with the 
normal (gaussian) distribution described by the moments Solti): Sq, (t1) 
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and S5,(t,) has to be used here. Then one repeats the procedure like in the 
point (2). 

4. The final distribution f(q,p,t,) at the time tn = tę + nAt is given by 
the end positions of the N randomly distributed points (q%, p% ). 

The statistical distribution will be more and more smooth with the 
growing number N of random points (q”, p”). 

When the time interval At between the sequential toss is small enough 
one can solve the set of eqs (8-12) approximately. Namely, it is possible to 
assume that locally around the point q the potential U is parabolic and the 
transport coefficients m,y and D are constant. These approximations allow 
to find the analytical solutions of eqs (8-9) and eqs (10-12) which become 
uncoupled. The solutions of eqs (8-9) after the i-th iteration reads: 


v v F A -a B - 
Gin. = G + (3 — rzą 14t _ at a2 At) jm (15) 
Pia = Ae! 484 Be“? (16) 
where 
pa U 
dq 
we JA. 
m 
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The coupled system of eqs (10-12) for the second moments has now the 
following solution: 
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For the small time interval At the right hand sides of equations (15-19) can 
be expanded in the Taylor series: 


dą = + At+= i= - PAŁ + 


1 
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Sy = — (At - ZAŁ +...) (23) 
y? 


Sy, = 2D[At - dh # (22, = w2JAĆ + ...] (24) 
The terms of the order higher than 3 in At are neglected here. 

The Fokker-Planck equation solved by the Monte-Carlo method is 
essentially equivalent to the Langevin equation with the normally distributed 
random force [10], when one omits in the equations (20-24) the terms of the 
order higher than 1 in Aż. 
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3. RESULTS OF THE CALCULATIONS 


The calculation was performed for the hot, fast rotating '®’Ir nucleus, 
which was experimentally investigated in [11]. We have taken the following 
Ansatze for the transport parameters when solving the FPE. The fission 
barrier of 1®7Ir was evaluated with the temperature dependent liquid drop 
parameters estimated in [12] The relative distance between fission fragments 
was taken as the collective coordinate and the fission barrier was minimized 
versus the neck parameter. We have got e.g. the fission barrier of the height 
Ug = 1.41 MeV for the nucleus temperature T = 2.5 MeV and the angular 
momentum J = 55h. The collective inertia m(q) in eq. (1) was assumed to 
be equal to reduced mass for relative motion of fission fragments and the 
friction parameter y was taken constant and it is given by the equation 


7 ee 
=—s=l.l 4 
8 0°*s 


The diffusion parameter D is related to the friction function by the Einstein 
relation: D = 7T. 

It is not easy to compare the whole distribution functions f(q,p,t) 
obtained by different method of solving of the FPE (1). So we have chosen 
a global factor related to f(q, p,t) 


F(t) = i 5 dq dp f(q,p,t), (25) 


which gives the probability that the deformation of the nucleus is smaller 
than that at the scission point qse. The number of nuclei which do not fission 
up to the given time t is proportional to F,,(t): 


N(t) = No Fsc(t), (26) 


where No is the initial number of nuclei at the time t = 0. 

We have solved the FPE within the local harmonic approximation (eqs. 
15-19) and using the Monte-Carlo method. The 10000 random points were 
generated. The initial shape of the distribution function f(q,p,t) was taken 
in the form of the probability distribution of the harmonic oscillator [13] 
which approximates the potential U(q) and mass m(q) around equilibrium 
deformation. The probability Fse (eq. 25) that a nucleus does not fission 
is plotted in Fig. 2 as a function of time. The different curves in the plot 
correspond to various time intervals At; between the sequential toss of the 
random points. Here, we have taken At; = (0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 
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t (1077's) 


Fig. 2. The Monte-Carlo estimate of the probability ¥,, (eq. 25) that a nucleus do not 
fission, At; = (0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6)1077) s 


WITr [= 55h, T = 2.5MeV 


0 10 7 20 30 40 
| t(107?'s) 


Fig. 3. The same as in Fig. 2. but for the first order in At expansion when solving the 
FPE, At, = (0.03, 0.05, 0.1, 0.15, 0.2, 0.3, 0.4)-1072! s 
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0.6)-107?! s. As it is seen from the figure we have got the convergence for 
At; < 107? s. 

Similar Monte-Carlo results as in Fig. 2 but for the expansion of eqs. 
(15-19) up to the first order terms in time are plotted in Fig. 3. Note that 
in this case the FPE is equivalent to the corresponding Langevin equation. 
We have chosen the following time intervals At; = (0.03, 0.05, 0.1, 0.15, 
0.2, 0:3, 0.4)-1072! s. The convergence was reached for At; < 0.5 - 10722 
s. The results of the Monte-Carlo calculation of Fse are compared in Fig. 4 
with the exact solution (solid line) of the FPE obtained by the propagator 
method (eq. 13-14). The dashed line represents the results from Fig. 3 
while the dotted line that from the Fig. 2. The both Monte-Carlo curves 
were evaluated for the time interval AT; = 0.5 - 1072? s, while in the case of 
the propagator method At; was four times longer. 


„0 - - - : - 
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t (1077! s) 
Fig. 4. The results of the Monte-Carlo calculation of Fse compared with the exact solution 


(solid line) of the FPE obtained by the propagator method. The dashed line represents 
the results from Fig. 3 while the dotted line that from the Fig. 2 


We can conclude that the the results obtained with these three methods 
are close to each other and that the numerical effort by the propagator 
method is the largest one. Additionally, we have found that this effort 
grows significantly with the dimension of the collective coordinates space 
while it is not the case for the Monte-Carlo method based on the first order 
expansion in At. 
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STRESZCZENIE 


Przedstawiono sposób rozwiązania równania Fokkera-Plancka dla jednej współrzędnej 


i sprzężonego z nią momentu pędu. Rozwiązanie otrzymane metodą propagatorów 
porównano z rozwiązaniem uzyskanym metodą Monte-Carlo. Pokazano, że obie metody 
dają zbliżone rozwiązania. Wykazano również, że równanie Fokkera-Plancka jest równo- 
ważne równaniu Langevina. 


